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Abstract 

Hotwire anemometers are used to measure instantaneous velocity from which the mean velocity and 
the velocity fluctuation can be determined. Using a hotwire system, it is possible to deduce not only the 
velocity components and their fluctuation but to also analyze the energy spectra and from that the 
turbulence length scales. In this experiment, hotwire anemometry is used to measure the flow field 
turbulence for an array of film cooling holes. The objective of this paper is to document the procedure 
that is used to reduce the instantaneous velocity measurements to determine the turbulence length scales 
using data from the film-cooling experiments to illustrate the procedure. 

Introduction 

Hotwire systems have several advantages for turbulence research including lower cost than other 
systems such as laser Doppler anemometers, high frequency response, small size, and lower noise/signal 
ratio (Ref. 1). Hotwire anemometry relies on the principle of convective heat transfer to measure velocity. 
The heated wire or film element is placed in the flow field; changes in fluid velocity will change the 
resistance of the wire. The wire being very thin instantly responds to the change in fluid conditions. The 
output of the hotwire is a voltage signal that can be calibrated to velocity from which mean and 
turbulence quantities can be determined. 

The objective of the current study is to characterize the turbulence field for an array of film cooling 
holes at a high blowing ratio where the jet is detached. The detailed turbulence measurements are an 
enabler for the development and improvement of numerical methods and tools for predicting film cooling 
flow and heat transfer. Turbulence is a familiar notion but remains a complex physical phenomenon to 
model and predict. Bradshaw (Ref. 2) defines turbulence as “a three-dimensional time-dependent motion 
in which vortex stretching causes velocity fluctuations to spread to all wavelengths between a minimum 
determined by viscous forces and a maximum determined by the boundary conditions of the flow.” Some 
describe turbulence in terms of an energy cascade in which large scale vortices are dissipated into smaller 
vortices which in turn dissipate into yet smaller scales and so on until the viscous effects dominate and the 
turbulence energy is dissipated into heat. A detailed description of the associated length scales and 
mechanism of turbulent transport is beyond the scope of this technical memo. The authors refer readers to 
several seminal texts on the subject (Refs. 2 and 3). 
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Experimental Apparatus 

The apparatus, shown in Figure 1, consists of an open-loop wind tunnel with a temperature-controlled 
coolant loop. The tunnel is a suction type tunnel that draws air from the room and passes it through flow 
conditioning screens prior to entering the test section. The freestream velocity is 30 ft/sec, and the 
Reynolds number based on freestream velocity and coolant hole diameter is 1 1,000. The test section is a 
square section 8.2 by 8.2 in. and 34 in. in length. 

The model consists of a three-hole array of film cooling holes that are fed by three long tubes. The 
model is inserted into the test section of the tunnel such that the heat transfer surface forms the wind 
tunnel floor. The cooling holes have a diameter D = 0.75 in., spacing P/D = 3, and are angled 30° from 
the streamwise direction. The coolant was fed through a manifold into three separate flow meters, then 
through 18 in. of hose and 12 in. of acrylic tube. The entire coolant path was insulated. 

Velocity and turbulence are measured using a two-wire hotwire system for obtaining components. 

The hotwire is powered by a constant temperature anemometry system, and data is acquired at a rate of 
50 kHz for 3 sec. Two different X-wire probes are used; one is used to measure the U- and V-components 
of velocity, and the other is used to measure the U- and VT-componcnts. In both cases, the two wires of the 
probe are separated by some distance and crossed to make an “X” shape. Figure 2 shows the probe 
(TSI Model 1246) used to measure the U- and VT-components, and Figure 3 shows the probe (TSI Model 
1240) used to measure the U- and V-components. 



Figure 1. — Wind tunnel test apparatus (left) and film cooling model (right). 


Model 1246 Cross Row “X” Probe, Sensors Upstream 



38 mm (1.S0) 
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Figure 2. — TSI Probe Model 1246 for measuring U- and l/T-components of velocity. 
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Model 1 240 Standard Cross Flow “X” Probe 


12.7 mm (.50) Hoi Film 
9.5 mm (.38) Hot Wire 


38 mm (1.50) 



Figure 3. — TSI Probe Model 1240 for measuring U- and ^-components of velocity. 

Procedure 

Calibration 

The probes are calibrated in-situ, with the probe located upstream of the film cooling holes at mid- 
height of the tunnel and the coolant flow shut off. The freestream velocity, measured using total and static 
pressure probes, is then varied from 0 to 150 ft/s, and the wire voltage is measured and recorded at each 
freestream velocity. A calibration curve for one of the probes is shown in Figure 4. The points are curve 
fit using a 4 th order polynomial that relates the delta voltage (the difference between the mainstream 
velocity voltage and the zero velocity voltage) for each wire to the freestream tunnel velocity. This 
velocity is assumed to be in the streamwise direction with no V- or VT-componcnts. The calibration is 
verified at the start of each run and repeated as needed. 

Velocity Trace 

During tunnel operation, the wire voltages are recorded to a data acquisition system and converted to 
velocity using the calibration curves to obtain the velocity components. As an example, Figure 5 shows 
the velocity trace at a location upstream of the cooling hole leading edge along the centerline of the tunnel 
and 1.0 in. above the floor of the tunnel. The mean velocity, (/ mean , is 32 ft/s. The standard deviation, u, of 
the velocity is 0.32 ft/s. Therefore the turbulence intensity, «/(/ mcan , at this location is 1 percent. 

Turbulence Energy Density Spectra Analysis 

The velocity signal shown in Figure 5 is in the time domain. By using a fast Fourier transformation 
(FFT), the data can be transformed to the frequency domain. The Mathworks Inc. software package 
Matlab is used for the analysis. The Appendix contains a listing of the Matlab script, parts of which are 
repeated here for clarity. 

FFT Analysis Transforms U From Time to Frequency Domain 

The velocity trace in Figure 5 contains 150,000 points (sampling frequency of 50 kHz for 3 sec). Let 
U be the array consisting of the 150,000 velocity points sampled over the 3-sec duration. ciFData is the 
FFT of the (/-velocity (i.e., the transform into the frequency domain). L is the length of the signal. 

a F Data = fft((7) 

L = length (aFData) 

The FFT of the (/-velocity signal results in half the number of unique points, so let 

n = L/2 = 150,000/2 = 75,000 
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Figure 5. — (/velocity trace in the freestream. 
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Plot Amplitude and Power Spectra of [/-Velocity 

To plot the velocity in the frequency domain, the absolute value of the FFT is used to calculate the 
magnitude, or amplitude, of the velocity array, i.e. 

aFMag = abs {aFData{ 1 \n)iL) 

Flere, L is used to normalize the data. The corresponding frequency array is 

aFreq = Fs*(l:ri)/n 

Figure 6 is the amplitude spectrum of the [/-velocity component. The amplitude ranges from 0 to the 
mean velocity 32 ft/s. The first point in the FFT corresponds to the zero frequency and mean velocity. 
This first point has been removed from the plot to enable us to analyze the turbulence energy. 


Single-Sided Amplitude Spectrum of U(t) 
0.025 


0.02 - 


0.015 - 



Frequency (Hz) 


Figure 6. — Amplitude spectrum of [/-velocity (omitting mean velocity at 0 Hz). 
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Power is defined as: 


Power = abs(fft(L0)- A 2/L 

Dividing by L ensures that the power is independent of the length of the signal. The energy E(f) can be 
determined by dividing power by the frequency. The energy density spectrum is shown in Figure 7. 

Integral Length Scale 

The integral length scale, A x (macro-scale), is a measure of the largest eddy size in the turbulent fluid 
and is calculated using the method shown by Roach (Ref. 4): 

[£■(/) [/mean 


where u is the standard deviation of the (/-velocity. To get E(f) as the frequency approaches zero from the 
data, the first 100 points are averaged, where the data asymptotes to a fixed value at the low frequencies. 
The dashed yellow line in Figure 7 shows E(f) as /approaches zero, and is determined to be 1.3445e 4 . 
Thus the length scale is calculated to be 

(1.3445e -4 ^— )(31.8 

A x = — = 0.0105 ft = 0.1256 in. 

4 (°- 322 §r) 

The length scale at this location (freestream inlet) is found to be about l/8 th in., which is roughly the size 
of the wire mesh upstream of the test section. 

Energy Density Spectrum 

10 ’ ^ 



Frequency (Hz) 


Figure 7. — Energy density spectrum of (/-velocity. 
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Dissipation Length Scale 

The dissipation length scale, X x (micro-scale), is a measure of the average dimension of the eddies 
responsible for the dissipation of turbulent energy. It is calculated as shown by Roach (Ref. 4): 


4 


2n 2 

U 2 u 2 


p CO 

f 2 E(f)df 

-'o 


The trapezoid rule is used for the numerical integration; the integration yields 2.3067e 6 ft 2 /s 4 . 


2n 2 


) 2 ft 2 ft 2 

x (31.82^X0.322^-) 


2.3067e' 


,A 2 


= 4.4001 e 5 — 

ft 2 


=0.0015 ft = 0.0181 in. 


Results 

The same procedure for calculating the length scale is repeated at each measurement location. 

Figure 8 shows the integral length scale along the centerline of the central cooling hole. Figure 9 shows 
the dissipation length scale along the jet centerline. The integral length scale is the size of the largest 
eddy. Note that at the hole location, the integral length scale is roughly 0.75 in. which is the diameter of 
the coolant hole. The dissipation length scale is smaller and is independent of the geometry. Whereas the 
integral length scale decreases in the flow direction as the coolant mixes and diffuses in the freestream, 
the micro-scale does not decrease or increase with the flow direction. The macro-scale decreases as the 
larger eddies break down into smaller eddies. 


4 


MacroLS [inch]: 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 



x [inch] 


Figure 8. — Integral length scale contour along jet centerline. 
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Figure 9. — Dissipation length scale contour along jet centerline. 


Conclusions 

The procedure for calculating the turbulence length scales from using hotwire anemometry is 
outlined, and a data point in the freestream inlet is calculated to illustrate the process. The instantaneous 
velocity signal obtained using a hotwire anemometer is processed using a fast Fourier transform to 
convert the signal from the time domain to the frequency domain. The magnitude is plotted against the 
frequency to yield the amplitude spectrum. The energy spectrum is obtained and plotted against 
frequency. It is from the energy spectrum and equations derived in Roach (Ref. 4) that the integral length 
scale (macro-scale) and dissipation length scale (micro-scale) are calculated. The results for the length 
scales along the centerline are plotted and discussed. The integral length scale is commensurate with the 
geometry. At the inlet, the integral length scale is approximately the size of the wire mesh upstream, 
while at the hole discharge, the integral length scale is roughly equal to the hole diameter. The integral 
length scale decreases in the flow direction as large eddies dissipate into smaller eddies. The micro-scale 
is at least one order of magnitude smaller than the integral length scale, and unlike the macro-scale, it 
appears to be independent of geometry and flow direction. 
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Appendix — MATLAB Script 


% Read in the two wire signal and process it to 


El = data ( : , 1 ) ; 

%let 

column 1 

be voltage 

E2 = data ( : , 2 ) ; 

%let 

column 2 

be voltage 

E0_1 = 4.667; 

%Zero 

voltage 

for wire 1 

E0_2 = 4.336; 
delEl = El - E0_1; 
delE2 = E2 - E0_2; 
al=0 . 1839; 
bl=-0 .198; 
cl=3. 9795; 
dl=l .8669; 
a2=0 . 134; 
b2=0 .259; 

%Zero 

voltage 

for wire 2 


get velocity 

1 (wire 1) 

2 (wire 2) 


c2=2 . 9707; 
d2=l . 8816; 

vell=al *delEl . A 4+bl*delEl . A 3+cl*delEl . A 2+dl*delEl; 
vel2=a2*delE2 . A 4+b2*delE2 . A 3+c2*delE2 . A 2+d2*delE2; 
U=0 . 5* (vell+vel2 ) ; 

W=0 . 5* (vell-vel2 ) ; 


Plot the U-velocity trace 


Fs=50000 ; 

L=length (U) ; 
t= (0 : L-l ) / (Fs) ; 
plot (t, U) ; 
title (' U-velocity trace'); 
ylabel ( 'U-velocity (ft/s) ' ) ; 
xlabel ( 'time (s) ' ) ; 


Sampling frequency of 50 kHz 
Number of samples 
Time vector 

Plot the velocity trace 


% Perform FFT analysis on the U-velocity signal and plot Amplitude spectrr 

aFData=fft (U) ; % Take FFT of U signal 

n=L/2; % FFT will yield half the number of unique points 

aFreq=Fs* ( 1 : n) /n; % Frequency array (half the length of signal) 

aFMag=abs (aFData (1 : n) /L) ; % Normalized Magnitude array (half the length c 

signal ) 
figure ; 

semilogx (aFreq (2 : n) , aFMag (2 : n) ) % Plot frequency against magnitude 
title (' Single-Sided Amplitude Spectrum of U(t)') 
xlabel (' Frequency (Hz)') 
ylabel ( ' I U ( f ) | (ft/s) ' ) 


% Plot the energy density spectrum 

Power=abs ( f ft (U) ) . A 2/L; %Power is the magnitude squared by L 

Energy=Power/Fs ; 

figure 

loglog (aFreq, Energy (2 : L/2+1 ) ) 
title (' Energy Density Spectrum') 
xlabel ('Frequency (Hz) ') 

ylabel ('Energy, E(f) (=Power/f requency) (ft2/s)') 
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% Calculate the integral length scale (macro scale) Roach Eqn 15 
uprime=std (U) 
uprimesquared=uprime A 2 ; 

Ubar=mean (U) 

Tu=uprime/Ubar 
Ef 0=mean (Energy (2:100) ) 

IntegralLengthScale=Ef 0 *mean (U) / ( 4 *uprimesquared) 

% Calculate the dissipation length scale (micro scale) Roach Eqn 8 
for i=l : length (aFreq) 

Y (i) =aFreq (i) * aFreq (i) * Energy (i) ; 

end 

Z=trapz (aFreq, Y) % Use trapezoid rule for integrat 

Z=Z*2*pi A 2/ (Ubar A 2*uprimesquared) ; 

DissipationScale=sqrt (1/Z) 
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